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Distinguishability of particles and its implications 
for peculiar mass transport in inhomogeneous media 
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A mass transfer directed from low to high density region in an inhomogeneous medium is modeled 
as a limiting case of a two-component lattice gas with excluded volume constraint and one of the 
components fixed. In the long-wavelength approximation, density relaxation of mobile particles 
is governed by a diffusion process and interaction with a medium inhomogeneity represented by 
a static component distribution. It is shown that density relaxation can be locally accompanied 
by a density distribution compression. In quasi one-dimensional case, the compression dynamics 
manifests itself in a hopping-like motion of diffusing substance packet front position due to a staged 
passing through inhomogeneity barriers and leads to a fragmentation of a packet and retardation 
of its spreading. A root-mean-square displacement reflects only an averaged packet front dynamics 
and becomes unappropriate as a transport characteristic in this regime. In a stationary case mass 
transport throughout a whole system may be directed from a boundary with low concentration 
towards a boundary with that of high one. Implications of the excluded volume constraint and 
particles distinguishability for these effects are discussed. 
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I. INTRODUCTION 

Mass transfer in inhomogeneous media and complex 
systems is often governed by nonclassical diffusion laws 
and exhibits peculiar transport effects. Specific examples 
include sub- or superdiffusion processes of, for example, 
Brownian particles in hydrodynamic flow [lJtl , and abso- 
lute negative mobility effect [a-[l3| . Very often, studying 
such effects we encounter a situation where relaxation 
goes in counterintuitive way, namely, instead of expected 
spreading density distribution undergoes compression, at 
least over certain time period, or even full collapse. In 
other words, mass transport is directed from low to high 
density region. In particular, such behavior occurs in 
systems with absolute negative mobility or negative dif- 
fusion coefficient. 

The phenomenon of absolute negative mobility or neg- 
ative diffusion coefficient may be of different origin. In 
particular, it can occur due to time- or particle-particle 
correlations, or nonequilibrium character of a process. 
This effect was demonstrated, for example, in electron- 
hole plasma, as an electron-drag effect arising from 
an electron- hole scatte ring J5j, in models of interacting 
Brownian particles |lll - fl3| or interacting lattice gas [lj] , 
and for nonequilibrium dynamics of Brownian particle 
in a periodic one-dimensional potential or a quasi one- 
dimensional channel in a presence of external force pe- 
riodic in time [7-9] . Note that the negative sign of dif- 
fusion coefficient entails absolute instability in a system 
and may signify the onset of a new phase formation, e.g., 
nuclei growth at first order phase transition [lj, [l5[ or 
collapse of electron-hole plasma that may be a precursor 
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for a neutral exciton formation, [lfj . 

The phenomenon of mass transfer from low to high 
density region also may be caused by an interaction of 
particles with medium inhomogeneity. In this case den- 
sity relaxation through this interaction may go much 
faster than via diffusion process, and, as a result, the 
latter may not be a main relaxation and entropy pro- 
duction mechanism [16]. A simple, although speculative, 
example of such a density relaxation is a non- wettability 
effect, caused by interaction between a liquid drop and 
a substrate (medium). Analogous effect in diffusive sys- 
tems will be of primary interest in this paper. 

Here we demonstrate that such an effect may naturally 
appear for diffusing particles in inhomogeneous media. 
In this case density relaxation dynamics in a subdiffu- 
sive regime can be locally accompanied by compression 
of density distribution. Moreover, in a stationary case, 
mass transport through the inhomogeneous sample can 
go from a boundary with low concentration to a bound- 
ary with that of high one. 

In order to demonstrate these effects we resort to the 
simplest model of a two-component lattice gas assuming 
that every lattice site can be occupied by one particle 
only (excluded volume constraint). In the limiting case 
of a frozen or static component this model reduces to the 
mass transport of mobile component in the inhomoge- 
neous medium. Another reason to use such a model is an 
appearance of a drag effects, which has been found in this 
framework in a presence of driving field and nonzero par- 
ticle hopping rates of both components, see, e.g., (l7l,ll8|. 
In addition, absolute negative mobility effect may appear 
in a case of non-Markov dynamics [10( . Two component 
lattice gas may also display memory effects. For exam- 
ple, correlation between jumps of tagged particle and its 
influence on tracer diffusion was examined in [19|, [20| . 
This correlation represents the tendency of a tagged par- 



tide to return to its previous position and it follows from: 
(i) a presence of different particle sort, i.e. particle dis- 
tinguishability, and (ii) the excluded volume constraint. 
The latter corresponds to infinite repulsion energy be- 
tween particles at the same lattice site. Formally, such 
"back correlations" describe a repulsion of a tagged atom 
from the other gas component. 



II. 



THE MODEL 



We consider the simplest model of a two-component 
lattice gas on a regular lattice, see, e.g., [19| . It is as- 
sumed that every lattice site can be occupied by a single 
particle of a sort m or n only. Particles can hop between 
neighboring sites separated by a distance a with rates v m 
and v n , respectively. Hopping events are assumed to be 
instantaneous, i.e., particle spends most of the time being 
localized on a lattice site. Rate equations for occupation 
numbers can be written as 119] 



drhi 
dhi 

~d7 



"mY^l^'ji 1 



rrii 



■ hi) — rhi(l — rhj — hj)], (la) 



v n 2_Jfij{l -hi- fhi) -hi(l- hj - rhj)]. (lb) 



Here m, = {0, 1} and hi = {0, 1} are local occupation 
numbers of m and n particles' sort at ith site, respec- 
tively, and sum is taken over nearest neighbors of ith site 
only. Equations ([T]) are based on the assumption that 
it is possible to neglect fluctuations in the number of 
jumps between sites i and j (2l|, and terms of the type 



j7Bj(1 



hi)dr gives a mean number of jumps (of 



m-particlc from site j to site i per time dr in this partic- 
ular case). 

In what follows we shall restrict our study to a macro- 
scopic dynamics. The simplest, if somewhat rough, way 
to obtain evolution equations for the average local occu- 
pation numbers rrii = (rhi) ((...) being statistical average) 
is to apply the mean field approximation 1221 or the local 



equilibrium approximation [2l|. The latter corresponds 



to the introduction of a local chemical potential associ- 
ated with a given lattice site, or a coarsened Zubarev 
statistical operator 23] . In this approach we lose all the 
information on fast (as compared to local equilibration 
time) processes and neglect any correlations. 

Next, we apply the long- wavelength approximation as- 
suming lattice constant a to be much smaller than char- 
acteristic inhomogeneity length scales l m and l n of m and 
n components, respectively, l m ~ l n 3> a. In the contin- 
uum limit, equations of motion for the densities take the 
form 



1 dm _,_ 

— = V Vm- 

v m dr 



{riS/m — mVn)], 



1 dn 
v n dr 



V[Vn — (mVn — nVm)], 



(2a) 



(2b) 



where we have introduced dimensionless coordinate r/a. 
Obtained equations describe smooth density profiles and 
can not be applied for length scales which are compa- 
rable with a lattice constant, l m ~ l n ~ a, where short- 
range correlations become significant, as near percolation 
threshold, see [19J ■ Note that in the case where inhomo- 
geneity length scales of two components drastically dif- 
fer, say l m 3> l n 3> a, higher orders of space derivatives 
should be taken into account in Eqs. ([2]). 

More general form of Eqs. © including terms taking 
into account external driving field have been obtained 
using mean field approximation |22J or phenomenological 
approach \YA_ 24 Ljpd exploited for investigation of phase 
transitions [13, [22, [24J, drifting spatial structures [25j | . 
and unusual transport effects in two-component driven 
diffusive systems [JM, (25[ . 

The main difference of Eqs. ([2]) from the ordinary 
diffusion equation is a presence of mixing flux j mn = 
nWm — m Vn = —] n m- This flux describes mutual drag 
of particles of one sort by particles of another one. It is 
caused by the two same reasons as a "back correlations" 
[l9j , mentioned above. The first one is a distinguisha- 
bility of two different sorts of particles, e.g., by spin or 
color. The second one is a local interaction (repulsion) 
between particles created by excluded volume constraint. 

The flows of gas components j m = —v m \Jm + 
v m (n\7m — mVn) and j„ = —v n \Jn + v n (mVn — uVm) 
can be represented as sums of diffusive and hydrodynamic 
parts. For example, for m-component, j, 



flow j^ 



j(l — n)Vm 



Jm+Jm>the 

-D m Vm describes diffu- 



sion of m-particles through vacant sites, unoccupied by n- 
particles, with local diffusion coefficient D m = v m (l — n). 
Term jj^ = m(—v m Vn) = mV m may be associated with 
transfer of m-particles by some hydrodynamic flow with 
velocity V m = -v m Vn, where concentration n of an- 
other component plays a role of a velocity potential. j„ 
can be analyzed along the same line. 

The drag of particles of one sort by particles of another 
sort directly follows from Eqs. ©. Consider the particu- 
lar case where characteristic inhomogeneity length scale 
of one component is much smaller than that of another 
component, say, l m 3> l n ^> a. To a first approximation, 
currents of both components are governed by a concen- 
tration gradient of the single component, Vn, at least for 
a certain period of time: 



1 dm 
v m dr 



V(m,Vn), 



1 dn 
v n dr 



V[(l-m)Vn]. (3) 



The drag effect particularly means that mass transfer 
of m-particles can be directed along their concentration 
gradient, Vm, i.e., from region with low concentration to 
region with that of high one. 

We will be interested in consequences of the drag ef- 
fect in the limiting case where hoping rate of n-particlcs 
is negligibly small in comparison with m-ones, v n <c v m . 
It is intuitively clear that this is some approximation of 
mass transport in an inhomogeneous medium where dif- 
fusion of mobile m-particles occurs through vacant sites 



unoccupied by "heavy" n-particles. Considering compo- 
nent n as a static one, we can write reduced equation of 
motion for a density of mobile component m 



dm 
~dt 



V[(l-n)Vm + TOVn] = (1 -n)V 2 m + mV 2 n, (4) 



where t = v m T. Equation (UJ) describes advection- 
diffusion with compressible flow (V 2 n ^ 0). Equation 
of this type is often used for the description of a trans- 
port in various systems [ll 0]. In a compressible flow, 
transport may exhibit intriguing effects as it is in the 
presence of stable foci (traps) of the flow |l| . In our case, 
the hydrodynamic flow, ~ mVn, is associated with a re- 
pulsion of mobile particles from the frozen component. 
Thus, density relaxation of mobile particles is governed 
by both diffusive mechanism and their interaction with 
"internal" field Vn. 



III. INFLUENCE OF AN INTERNAL MEDIUM 
FIELD 



which means that mass transfer may occur towards 
higher concentration region. Contrary to that, mobile 
particles tend to evacuate from regions close to maxima, 
where n(r max + Sr) ps n max (l — q 2 (Sr) 2 ), which can be in- 
terpreted as forcing mobile component out by the frozen 



m(r, t) w m(r, 0) exp(-tq 2 n max ). 



(7) 



Such exponential dynamics indicates a presence of the 
faster process than that of diffusive, which locally, over 
short time periods, leads to the dependence Ri oc ~ t, 
while in the case of anomalous diffusion R ~ i£ , with 
C < 1- In this regard, such global integral characteristic 
as root-mean-square displacement 



R{t) = vV) 



J x 2 m(x,t) dx\ 
J m(x, t) dx 



(8) 



may be improper for the description of a diffusion process 
due to loss of information on fast local dynamics. 



As it is well known, the advection-diffusion can exhibit 
anomalous behavior that manifests itself, in particular, 
in non-classical dependence of a root-mean-square dis- 
placement on time, R = (r 2 ) 1 / 2 ~ t^ , where £ ^ 1/2 
is the exponent of anomalous diffusion [1|. In the case 
of compressible flow considered here (V 2 n ^ 0), the ex- 
ponent C may lose universality and may depend on the 
degree of compressibility. Such anomalous phenomenon 
is not necessarily takes place asymptotically in the limit 
t — > oo, but can occur, at least, on a time scale that is 
less than a finite mixing time |l|, |j| . On this time scale, 
density relaxation may strongly depend on initial state 
and demonstrate peculiar behavior. 

In this section we show that mass transport from low 
density region towards a dense one in an inhomogeneous 
medium may locally accompany subdiffusive (£ < 1/2) 
process. Moreover, such integral characteristic as a root- 
mean-square displacement R does not describe relaxation 
process properly. Unlike in the case of ordinary diffusion, 
motion of packet front rj(t) does not coincide with a 
root-mean-square displacement time dependence and has 
phased behavior. 

Indeed, at relatively low values of gradient Vm mass 
transfer is determined by the second term in right hand 
side of Eq. ((4]), rh ~ mV 2 n, which, at least over short 
time period, leads to the dependence 



m(r,t) ~ m(r, 0) cxp(iV n(r)). 



(5) 



Behavior of a density profile mir, t) in the vicinity of 
frozen component distribution n(r) local maxima and 
minima is different. Near minima, where n(r m i n + Sr) sa 
n m i n (l + q 2 (Sr) 2 ) 1 mobile particles tend to accumulate, 
i.e., initial density profile is squeezing, 



A. Packet fragmentation and hoping dynamics of 
its front position 

In order to illustrate features of the transport phe- 
nomenon in an inhomogeneous medium described by 
Eq. Q, we consider relaxation process in quasi one- 
dimensional case, supposing that transverse size L±_ of 
a system is of the order of magnitude or less than char- 
acteristic inhomogeneity length scale of the densities m 
and n, L± < l m ~ l n . Then equation |4]) reduces to 



rh = (1 — n)d x m + md x n, 



(9) 



where d x labels one-dimensional derivative. 

We will consider spreading of initial Gaussian distri- 
bution 



m(x,0) = Mcxp(-x 2 /4l 2 ) 
in a periodic field of the frozen component 
n{x) = (JV/2)(l-cosfc a;). 



(10) 



(11) 



m(r, t) « m(r, 0) exp(ig 2 n r , 



(6) 



Degree of medium inhomogeneity is determined by the 
amplitude N and period 2irkQ of the frozen component. 
As it can be seen from Fig. QJa), presence of inhomo- 
geneity slows down packet spreading. Such slowing down 
may be caused by decreasing of effective diffusion coeffi- 
cient (renormalization of time) and/or decreasing of in- 
dex £, which characterizes subdiffusive regime (see inset 
at Fig. [Ha)). However, similar behavior of index £ (see 
inset at Fig. [TJb)) is given by Eq. ([9]) with only diffusion 
term included rh — d x [(l — n)d x m], where relaxation has 
another character, Fig. [D(b). Packet fragmentation is as- 
sociated with the second term in Eq. ((4|), which describes 
repulsion of light atoms from heavy ones. If this term is 
dominant, mass transfer does not go through a diffusive 
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FIG. 1: a) Density distribution m(x,t) for different values of 
N at t — 31888; b) mi(x) is the density distribution, given by 
equation mi = d x [(l — n)d x mi], rri2(x) obeys Eq. ((9]). Insets 
display corresponding index behavior C(lni) = ln_R(t)/lnt, 
where R(t) is given by Eq. ([8]). Initial amplitude M — 0.1, 
initial half- width \f% — 14.14, and fco =0.1 are used in both 
(a) and (b) cases. 



mechanism, described by the first term in right hand side 
of Eq. (T3|) , but instead occurs due to interaction with a 
heavy subsystem. 

As it directly follows from Eq. (TJ|, local flow j m of the 
mobile component is defined by two contributions, 



-(1 — n)d x m — md x n. 



(12) 



The first term, j^ = —(1 — n)d x m ~ —(1 — n)(m/l m ), 
where l m is typical length scale of m{x, t) variation near 
x, describes standard diffuse spreading, i.e., tendency of 
transfer to be directed towards concentration lowering. 
Presence of the frozen component leads to a decrease 
of diffusion coefficient by the factor of 1 — n. Frozen 
component plays the role of a barrier for mobile par- 
ticles. As concentration n(x) (barrier height) increases 
the velocity of penetration through the barrier by mo- 
bile atoms (diffusive flux j^J decrease as 1 — n. The flux 
j^j = —md x n ~ —m(nko) describes repulsion of mobile 
atoms from frozen one. Fluxes j^ and jj^, defined by a 



density gradients d x m and d x n, may have different signs 
and behave as two competing fluxes. If 1/n — 1 < l m ko, 
then total flux j m may be locally directed towards region 
with higher concentration, which means local compres- 
sion of mobile atom distribution. 

Figure [5] illustrates dynamics of initial Gaussian dis- 
tribution spreading, which has phased nature. Front 
propagation through the first barrier towards neighbor- 
ing local profile minimum (depicted at Fig. H£a)) con- 
sists of two stages. At first, particles locally accumulate 
in the minimum due to jj^ domination until the condi- 
tion 2/N — 1 > l m ko for m(x,t) profile is met, or in 
other words, until local "pressure" (ex d x m) gets high 
enough to overcome repulsion from next barrier. Expres- 
sion 2/N — 1 > l m ko defines the condition for diffusive 
penetration through the barrier. Then process repeats 
for the next local minimum, while in the previous one 
particles undergo usual spreading, Fig.^b). Such staged 
process describes motion of a packet front Xf(t), which 
position is determined by condition m(xf,t) — 0, see 

Fig.rjTc). 

Front Xf(t) divides x-axis in two regions: (i) \x\ < Xf, 
where concentration m(x,t) is decreasing for any x and 
density relaxation is going mainly due to diffusion with 
diffusion coefficient renormalized by a profile of the frozen 
component; (ii) \x\ > Xf, where concentration is increas- 
ing and relaxation is strongly affected by a repulsion from 
the frozen component (medium inhomogencity) . The lat- 
ter leads to a local compression and density increase. 

As can be seen from Fig. HJc) , medium inhomogeneity, 
which in our case is determined by an amplitude N , leads 
to appearing of difference between packet's front motion 
Xf(t) and root- mean-square displacement R(t) dynamics. 
Front Xf (t) exhibits phased behavior and is governed by 
two linear in time, fast and slow processes. Packet front 
spends most of the time being localized inside a barrier 
and performs quick jumps into neighboring barriers while 
root-mean-square displacement (Eq. ([5])) gives only an 
averaged dynamics of front motion, see Fig. [He) . 



B. Compressibility of initial density distribution 

Compression (growth) of density distribution also can 
occur for the central packet during certain initial time 
interval. Indeed, from Eqs. (|4")) . (ITU)) . and (TTTJ) it is easy to 
estimate behavior of distribution amplitude at the peak 
over short time periods, 



m(0,t) « M + ^(Nk 2 l 2 - 1) 



(13) 



Amplitude m(0,t) tends to increase if N > (kol)~ 2 . On 
the other hand, density distribution compression means 
that velocity (time derivative) of a mean-square displace- 
ment v(t) = d(x 2 )/dt is negative. The latter may be 
roughly estimated in Fourier domain 
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Packet dynamics during time intervals: 



[to,ti 



FIG. 

b) [ti,t 2 ]; (to = 0, ti = 14 X lO 13 , t 2 = 45 x lO 13 time units). 
Initial amplitude M = 0.5, amplitude of the frozen component 
N — 0.7. Shaded pattern shows density distribution profile of 
the frozen component n(x), given by Eq. (|lll) . arrows denote 
whether sub-packet is squeezing (growing) (f) or spreading 
(4) on certain time interval. Dashed lines Xf denote region of 
main localization during corresponding time interval. 

c) Time dependence of packet front Xf(t) defined as a point 
where m = 0, and a root-mean-square displacement R(t) de- 
fined by Eq. 0, for different values of N. Lines' labels cor- 
respond to values of N. Dotted lines x^ correspond to the 
dotted lines at Figs. [2] (a) and (b). 




FIG. 3: Central peak squeezing (growth) of a mobile com- 
ponent density profile during initial time interval [£o,£i], 
(to = 0, (i = 3x 10 4 time units). Initial amplitude M = 0.1, 
N = 0.9, initial half-width \/2l = 1.7, and k = 1.05. Shaded 
pattern shows density distribution profile of the frozen com- 
ponent n(x), given by Eq. (|lip . arrows indicate directions of 
packet evolution. Filled region at the inset is determined by 
inequality (|19[) and condition N < 1, and corresponds to a 
negative velocity of packet motion. 



where m k (t) is the Fourier image of m(x, t) and obeys 
equation 

m fc = -(N/A)[(k 2 - 2kk )m k _ ko ] - 
-(N/4)[(k 2 + 2kk )m k+ko ] -k 2 [l- (N/2)]m k . (15) 

Since v(t) is determined by limit fc — >■ 0, it is reasonable 
to approximate Eq. (|15[) by its expansion for small k (k <C 
feo) 

m k » -k 2 [l - (N/2)]m k - (N/2)k 2 m ko 
~Nk 2 k a —m ko - (l/A)Nk 4 -^m ko , (16) 

where we have omitted terms with vtii ka ~ assuming 
that typical scale of inhomogeneity for m(x, t) is of the 
order or less than that of for n(x), i.e., k < kg. Taking 
limit k —y we obtain expression for velocity 



vns(2-N)-—(l + 2fc J- )m ko . 
m \ dk / 



Here m k „ is defined by the equation 



m ko 



N 



N 



M 1_ t m *o + ^r m o 



(17) 



(18) 



From (|17p condition for packet squeezing (v < 0) directly 
follows, see inset at Fig. [3J 



2;2\ 



2/N < 1 + (4fcgr - 1) exp(-/sgr) 



(19) 



Numerical solution of Eq. ^ corresponding to such a 
case is shown at Figure [3] Initial compression lasts for 



certain time period after which diffusion goes by means 
of sequential overcoming of potential barriers (as it was 
discussed in previous subsection). Note that the estimate 
(j!9[) is rough and initial packet squeezing actually appear 
at smaller values of amplitude N, as shown at Fig. [3] 



V. 



SUMMARY AND DISCUSSION 



IV. 



THE STATIONARY CASE 



As it was shown in the previous section, interaction 
of diffusing particles with inhomogeneities of a medium 
leads to a local particle accumulation, i.e., local mass 
transport towards higher concentration. It is interesting 
that such transport effect can be realized at global scale, 
throughout a whole system, with arbitrary distribution 
n(x). In order to illustrate such a possibility we consider 
quasi one-dimensional boundary problem 



(I — n)d x m + md x n = 



with boundary conditions 



to(0) = m(0), m(L) = m(L), 



(20) 



(21) 



where L is a sample length. Solution of Eq. (I2U1) has the 
form 



, m(0) 

mix) = (1 — n(x))\ - 7^7 — J 



< 



l-n(0) J (l-n(0) 



. (22) 



where J is a total particle flux through the system 



J 



<% 



o [l-^)] 2 



m{L) 



m(0) 



l-n(L) l-n(0) 



= -J / jm(») dx. 

L Jo 



(23) 



In the last expression j m (x) is given by Eq. ([12]). 

As we see, stationary density distribution m(x) may be 
inhomogeneous even if concentrations at the boundaries 
rn(L) = m(0) are equal. 

As it was mentioned above, mass transfer in such a 
system is governed not only by the mean field of m(x) 
gradient, but also by interaction of a mobile subsystem 
with a frozen one. As it follows directly from Eq. (fT2|) 
mass transfer may go from lower concentrations towards 
higher one. Assuming m(L) > m(0) one gets, that flux J 
is directed towards the higher concentration region, i.e., 
from the boundary x = to the boundary x = L, if the 
following condition is satisfied 



m{L) l-n(L) 
m(0) < 1 - n(0) ' 



(24) 



Thus total flux direction is determined by boundary con- 
ditions for m and n components. 



The presence of a second sort of particles in lattice gas 
leads to peculiar transport effects. In the case of a two- 
component gas mass transport is affected by an action of 
an additional flux, contrary to the case of a single com- 
ponent gas of indistinguishable particles. This flux is as- 
sociated with mixing of different gas components (inter- 
diffusion term) and in a presence of interaction between 
particles such as excluded volume constraint in our case, 
and may lead to a drag of one sort of particles by another 
one. If particles of one sort have negligibly small mobil- 
ity in comparison with another sort, so that one sort is 
assumed to be static, then the mixing flux, or its part, 
of mobile particles transforms into a stationary hydrody- 
namic flow that may drag mobile particles. As a result, 
density relaxation of mobile particles occurs via diffusion 
process and/or via an interaction with a medium, i.e., 
with a frozen gas component. 

In our paper we have presented a case of a system 
where diffusive mechanism cooperates with density re- 
laxation under action of "internal medium field" , given 
by a periodic distribution of a static component. Such 
mutual cooperation leads to unusual transport with local 
mass transfer resulting in an increase of concentration. 

Indeed, as it was shown, if medium is strongly inho- 
mogeneous, its "field" leads to a local accumulation of 
diffusing particles in minima of inhomogeneity profile. 
Presence of such a local processes is manifested by frag- 
mentation of a packet during its spreading. The dynam- 
ics of a diffusive process also becomes peculiar and mass 
transport occurs as sequential penetration through inho- 
mogeneity barriers, i.e. regions with high concentration 
of the static component. This leads to a phased character 
of a packet front motion that consists of a sequence of fast 
and slow processes (moves) with linear in time coordinate 
dependence. Such front behavior is similar to motion 
of a defect (slow motion inside barrier and fast passing 
through inhomogeneity minima). Contrary to the case 
of ordinary diffusion behavior of a root-mean-square dis- 
placement R(t) and that of a packet front Xfit) do not 
coincide. R(t) reflects only an averaged packet front dy- 
namics. Local accumulation of particles, which accom- 
panies diffusive process and sequential barriers overcom- 
ing, also slows down packet spreading. The latter entails 
subdiffusive regime. In stationary case interaction of dif- 
fusing particles with medium inhomogeneity can lead to 
mass transport directed from a boundary with low con- 
centration to a boundary with that of high one. 

Note, that quasi one-dimensional case considered in 
this work is rather illustrative, because dimensionality re- 
duction leads to enhancement of particle-particle correla- 
tions. In two- and three-dimensional cases demonstrated 
effects may not be so pronounced. However, accounting 
for interparticle interaction on nearest neighboring sites 
leads to essential changes in mass transport [l4j and may 
enhance discussed effects. 

Note, that assigning appropriate mobilities to compo- 



nents in Eqs. ([2]) we can model not only inhomogeneous 
media but also systems with quenched randomness, as it 
was pointed out in [lj, [24|, or, for example, quenched 
inhomogeneities, dynamically generated in glasses |26| . 

Distinguishability of two sorts of particles and ex- 
cluded volume constraint are responsible for unusual 
transport effects appearing in two-component lattice gas. 
Note, these effects may speculatively be considered as a 
purely statistical consequence of particles' distinguisha- 
bility in a multicomponcnt non-interacting Fermi gas. 

Equations ^ have been obtained using rough approx- 
imations, which entails certain limitations on their valid- 
ity. First of all we neglect fast processes in the system, as 
it is generally the case for the diffusion approximation, 
the framework basic Eqs. ([T]) are written in. In addition, 
Langevin source of fluctuations in a number of jumps 
between lattice sites was omitted in Eqs. ([TJ. We also 
loose information on fast processes applying the mean- 
field approximation which is equivalent to the local equi- 



librium one. This approximation also means that we ne- 
glect any correlations in the system, in particular, short- 
range "back correlations" (local memory effect) [19|, [2CJ 
that are known to contribute to effective diffusive pro- 
cess. The long- wavelength approximation does not allow 
to consider mass transport on length scales comparable 
with the lattice distance. In particular, it is in a vicinity 
of percolation threshold where short-range correlations 
become significant. 
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